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possible to formulate string theory nonperturbatively solely in terms of 
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§1. Introduction 

It is now widely recognized that the nonperturbative effect in string theory 
that behaves as ~ e^li^^^ stems from the dynamical hypersurfaces in space-time 
(the D-brane) on which open strings can end.^^ The double scaling limit of matrix 
models^)' ^) may enable us to learn more about this, at least for noncritical strings. It 
was realized in the early 1990s that the string equation (Painleve equation) indeed 
inherits such nonperturbative effects;^)' they correspond to the deviation of its 
solutions from the genus expansion of two-dimensional gravity. In the early 2000s, 
studies of the D-brane have been developed from the viewpoint of Liouville field 
theory.^)"^) 

Given this situation, we are tempted to reinvestigate nonperturbative effects of 
the noncritical string theory in full detail. In particular, we would like to ask if 
nonperturbative effects introduce a continuous parameter characterizing vacua like 
the 6 parameter in QCD, or if they are fully calculable as an intrinsic property of 
the matrix model. 

We focus on the c = noncritical string here in order to address the above 
questions as clearly as possible. First, let us recall the nonperturbative effects ob- 
served in string equation, following Refs.^^ and.^^ The variable for which to solve is 
the specific heat u{t), depending on the renormalized cosmological constant t. The 
string equation with respect to it takes the form of Painleve equation I, 

u^ + ^diu = t. (1.1) 

The leading-order contribution to u{t) in the genus expansion is (5s)°, and it is 
related to the free energy F {F = lnZ with the matrix model partition function Z) 
through 

u = gldlF. (1.2) 
u{t) admits a perturbative series expansion: 

»P.»(«) = -VI+j^| + _|_| + .... (1.3) 

However, there may be some degree of ambiguity in the above solution, due to 
nonperturbative effects. To examine this point, let us suppose that there is another 
solution u{t) infinitesimally close to the above perturbative solution, which we write 

u{t) 

— '^pert 

{t) + Au{t) . (1.4) 
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Then, if we write Au{t) in the form 

h{t) 

Au{t) = Ae~^, (1.5) 
the equation for h{t) admits an expansion with respect to gs, 

(^{dth{t)f - uVt + Ae-^^ + gs {-dfh{t)) + 0{g^s) = ■ (1.6) 

Noting that e is neghgible as long as gs is small, h{t) can be solved as 

^(i) = ^it + |lni + 0(5,2). (1.7) 

By inserting this into Eqs. (1.5) and (1.4) and integrating u{t) twice over t, the free 
energy F is found as 



— -^pert ~l~ -finst 

5 

1 n ^ 



^-- = -15^-48 ^^^+ 576o|+^(^^)' 

(l + 0(5.))- (1-8) 



C 

-Pinst = — exp 



8\/3 5 

■t4 
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Here, i^inst is the non-perturbative contribution to the free energy, which can be in- 
terpreted as an instanton effect.^)' Recent developments have identified the origin 
of these effects as D-branes.^*^^"^^-* However, the overall factor of C for the instanton 
contribution cannot be determined from the string equation (1.1) itself; it is a con- 
stant of integration that shoTild be determined from the other boTindary condition 
mentioned above. That is, the string equation determines the D-instanton action 

5inst = ^-^t\ (1.9) 

55s 

uniquely, and it shows that the D-instanton contributes to the free energy in the 
form 

— e"^'°=* . (1.10) 

However, the weight C relative to the perturbative contribution, which corresponds 
to the chemical potential needed to bring a D-instanton into the system, remains 
undetermined with the string equation (1.1). 

Our main focus here is concerned with the chemical potential of the D-instanton. 
There are the following two possibilities for C: 
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(a) (7 is a parameter charaeterizing the vacua as a result of quantizing the system, 
like the ^^-parameter in QCD. Specifically, each value of C corresponds to a distinct 
vacuum. 

(b) C is calculable and its value is uniquely determined. In this sense, the string 
equation does not fully describe the properties of the system. 

We are able to determine which possibility is actually realized by inquiring 
whether the chemical potential of the instanton can be calculated directly by carrying 
out the path integral of the c = matrix model. 

Our main result here is the following: The value of C can be calculated. More- 
over, C does not depend on the details of the matrix model action, and it is thus a 
universal quantity, given by 

C = i ^ . (1.11) 

8-33 • 

Reflecting the instability of the vacuum in the presence of a D-instanton, C is a 
purely imaginary number. Our result shows that the lifetime of such a vacuum 
is uniquely determined and does not depend on the regularization scheme. These 
results indicate that assertion (b) is correct. 

This paper is organized as follows. We start in §2 by identifying the instanton 
effect in the matrix model and showing that it corresponds to the D-instanton effect 
in c = Liouville field theory. Then, we divide the partition function of the matrix 
model into the sum of contributions ^("-^^t) from n of instantons {n = 0, 1, • • • ) 
by restricting the integration regions of the eigenvalues in a proper way. There, we 
also derive various significant equations for estimating the chemical potential fi of 
an instanton, i.e., the ratio of the single- instanton contribution Z^^~^^^^^ to the trivial 
vacuum contribution Z^^'^^^^\ 

^^•(l-inst) 

In §3, we apply the expression obtained in §2 to the estimation of the chemical 
potential of an instanton, with the help of the orthogonal polynomial method. We 
find that the chemical potential of an instanton is indeed given by the value in 
Eq. (1.11). 

Section 4 constitutes the second part of the paper. Then, we reconsider the result 
obtained in the previous sections from the matrix model calculation by examining the 
loop equations more closely. First, we consider the loop amplitude in the background 
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produced by singlc-instanton. This is done by taking the expectation value of the 
loop amplitude after restricting ourselves to single-instanton sector a priori, ignoring 
the other sectors. We find that the loop equations in the large N limit determine 
the loop amplitude in the single-instanton sector, up to a constant factor. In order 
to determine this constant, we need the loop equations at all orders. 

Next, we examine whether the loop equations determine the chemical potential 
of an instanton. As seen in §3, the chemical potential of an instanton is universal. 
However, it requires regularization of a divergence of type InA^. Thus, we find that 
the loop equations obtained after taking the continuum limit cannot determine the 
chemical potential of an instanton. This fact indicates that it is quite difficult in 
string theory with closed strings only to fully describe the nonperturbative efi"ects of 
strings, and that open strings or matrices should be regarded as fundamental degrees 
of freedom. Section 5 contains the conclusion of this paper. 

§2. The instanton in c = noncritical string theory 

In this section, we show that the instanton in the one-matrix model is identical 
to the D-instanton in the c = noncritical string theory and that it indeed gives the 
nonperturbative efi^ect discussed in the introduction. 
2.1. Action of the instanton 

As a concrete example, we consider the one-matrix model with a cubic potential, 

5 = A^try (</,), y(^) = 1^2 _ 1^3^ (2.1) 

where (f> is an N x N Hermitian matrix. The effective action for the eigenvalues Aj 
{i = 1, - ■ ■ ,N) is given by 

Ses = - ^log(Ai - Xjf + A^5^ F(Ai). (2.2) 

i<j i 

Hereafter, we consider the situation in which a single eigenvalue (Ajv) is separated 
from the others. Then, the partition function of the matrix model (2.1) is expressed 
as 

Zn = JdxJ dXi---dXN-i (jlix-Xif^ Z\(^^-^)(Ai,--- ,^N-if 
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where x denotes Ajv, and A^^^^\Xi, • • • , Xjsf^i)"^ is the Vandermonde determinant 
in terms of Ai, • • • , Aat-i. Thus we have ^flilli^l^ ~ -^j)^) A^^~^\Xi, ■ ■ ■ , Aat-i)^ = 
A^^\\i, ■ ■ ■ , Ajv)^- By introducing an (N — 1) x (A'' — 1) Hermitian matrix cp' , this 
can be rewritten as 

= Z'j,_, J dx <det(x - 0')")' e""^^^"^^ 

= Z'j^_i J dxe-^^-«^''\ (2.4) 

where 



{Oy = -y— / ^'Oe-^*^^^'^'). (2.5) 



In the large- AT hmit, we can set 



^ tr f -J— ) ) = R{x), (2.6) 



^N-l \x-(f)' 

where R{x) is the resolvent of the matrix model (2.1), which is given by 

R{x) = (^tv- ^ 



N X- 



= I + v'^'(x)2+p(x)) . (2.7) 

Here, p{x) is a polynomial of degree 1, and the branch of the square root is taken so 
that R{x) ~ 1/x as a: — oo. Therefore, in the large-iV limit, VeS becomes 



V}^\x) = -2Re f dx'R{x') + V{x) 

J X* 

= -Re / dx'^yV'{x'y +pix'), {21 

•J x» 



where fixes the origin of . Here, we note that when x' in the above integration 
is inside the region where the other eigenvalues are distributed, we have to take 
the real part of the resolvent, because the square of the characteristic polynomial 
det(a:' — 0')^ is originally positive-definite. Accordingly, {x) has a plateau in this 
region, as shown in Fig. 1. The meaning if this plateau is that the force from the 
other eigenvalues acting on the eigenvalue we are considering vanishes in this region. 
This is natural, because inside the sea of eigenvalues, the force should be balanced 
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to zero. Thus if we choose a value of x, x*, inside the plateau, V^g (x) gives the 
effective potential in the case that we regard its origin as the value at the plateau. 
In particular, its local maximum V^lf\xi^o)) can be regarded as the potential of the 
instanton that lies at the top of the effective potential V^q\ A simple computation 
shows that in the double scaling limit^ we have 

g = g^e^pi-eh), = exp(eC), N^ = e~^. (2.9) 

The height of the potential barrier then becomes 

^^(^^^(^(O)) - vS\x^,^)) = ^tt, (2.10) 

where t is normalized so that the 'specific heat' u derived from the free energy F(t) 
as u = F satisfies the Painleve equation, 

2 1 / N 

u^ + -u = t. (2.11) 



^ Here we have set the string coupling constant as = 1. thus, all dimensional quantities in 
the continuum limit become dimensionless in units of gs- In §2 and §3, we follow this convention. 
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Correspondingly, (" is normalized so that the disk amplitude obtained by taking the 
scaling limit of (2.7) becomes 

*(c) = (c-^vt) \f^+7t. (2.12) 

Note that (2.10) agrees with the prediction of the Painleve equation, eq:sinstpv).^) 
In §3, we examine how the coefficient as well as the power of t in (1.8) can be 
determined from the point of view of the matrix model. 
2.2. Chemical potential of the instanton 

Next we take a closer look at the contribution from the instanton to the free 
energy. We begin with the partition function Z]\f, which can be written as 

Zn = j dXi---dXNG, G = A^^\X^,---,XNfe^p(^-Nj2v{Xi)j. (2.13) 

This can be divided into the sum of the multi-instanton sectors in the large- AT limit, 
where the interactions between the instantons are of 0{1/N) compared to the leading 
contributions, as we see below. Thus we write 

Zn = 4°""'^ + 4'""'*^ + 4'""'*^ + • • • > (2-14) 

where the A;-instanton sector is characterized as that in which k eigenvalues are 
separated from others; that is, they do not lie in the region < x < X(2). 
Let us consider the 1-instanton sector: 

Zg-^t) =N dx H dXiG{x, Ai, • • • , A^_i) 

= NZ'j.J'-'"''''^ f dx {detix - <^')2)'(°-'"^*) e-^^(-) 

Jx<x^i'),X(^2)<3: 

= NZ'^J"-'^'^ [ dxfix), (2.15) 

Jx<X(^iy,X(^2)<3: 

where the overall factor N reflects the number of ways of specifying the isolated 
eigenvalue. Quantities with primes are defined as in (2.5): 

(0)'(°-'-*) = f dcp'Oe-^'^nfX (2.16) 
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Here, all eigenvalues of <p' are understood to lie between a;(i) and X(2), and 

f{x) = {det{x - e-^^W. (2.17) 

It is evident that in (2.15) if we change the interval of the integration with respect 
to X from {x < a;(i),X(2) < x} to {iC(i) < x < X(^2)}j this integral would give the 
0-instanton partition function multiplied by N. This observation leads us to the 
relation 

iVZ^_/°"°'*^ / dxf{x) = NZ^^'''''\ (2.18) 

Thus, we derive the chemical potential of the instanton in terms of the correlator of 
the matrix model (2.1) as 

(i-inst) J ^ dx(det(x-</))2)(™)e-^^(-) 

where we have used the fact that in the large- AT limit (or the double scaling limit) 
we have 

fix) = (det(x - </,)2)(o--«t) e-^^W, (2.20) 

as in (2.6). 

Similarly, "^^^^ is given as 

^(fc-inst) _ r< f 7' (0-mst)\'^ 
^JV — N<^k y^N-k J 

X n f / dx, (det(x, - ff)' ^'-'"'^'^ e-^^(-')) 

i=N-k+i yxi<^(i),x(2)<xi ) 

xZ\«(xAr-fc+i,--- ,^iv)', (2.21) 

where the prime indicates quantities of the (N—k) x (N—k) matrix model. As long as 
k <^ N, the last factor, Z\^^)(x7v-fc+i, • • • , a^Ar)^, can be ignored, because it gives only 
an 0{1/N) contribution compared to the leading one. Physically, this corresponds to 
switching off interactions between instantons (the dilute gas approximation). Then, 
the above equation becomes 

Zp-''^=r,cJz'^J'-'''''^y([ dxfix)] . (2.22) 

Repeating the same argument, we obtain 



nC^ iz'^-^'-'^''') ' I dxfix) = NC.zf'-'^^ 

\Jxn)<x<xi2) I 



(2.23) 
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and 

This implies that the free energy is in fact given by 

F = + //. (2.25) 

2.3. An instanton as a D-instanton 

In this subsection we confirm that the instanton of the matrix model we have 
considered is indeed the D-instanton in the c = noncritical string theory. 

We first note that in ordinary critical string theory, contributions from a D-brane 
correspond to adding surfaces with open boundaries. Let us check that this is also 
the case with our instanton. We rewrite the partition function (2.3) as 

Zn = JdxJ dc^'dqdq e-5(-'- )-W(x)^ (2.26) 

where 

^(i-inst) ^ ^trl^(c/>') + q0 - x)q,, (2.27) 

i=l,2 

and qi is Grassmann odd in the fundamental representation of U{N — 1). This 
implies that if we evaluate the integration over x by using the saddle point value at 
X = X(Q), the interactions over q and q provide new contributions from surfaces with 
open boundaries^^^"^^^ for the partition function. This can be expressed as 

' q' q q 

From these considerations, we are led to suspect that our instanton is identically the 
D-instanton. 

In order to make a precise comparison with the continuum (Liouville) theory, 
let us consider a loop amplitude in the instanton background. If we take account of 
the effect of the instanton, the resolvent can be written 

/I f I 

:tr 



N \z 
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1 



N-l 



\ 1=1 

N-l 

X 

i=l 



^(O-inst) , ^ / ^ N > (0-inst) ^(l-inst) / ^ ^ / , A^-l 

Zjv \iV \z-(t)) Zn \Nz-X(o) \N^z-Xi 

+ • • • , (2.29) 

where we have replaced the integration with respect to x over the interval x < 
X(2) < X with the saddle point value at x = a;(o), which is valid in the large- A?^ limit, 
and we have used 

' 1 ^"^ 1 \ 1 f / 1 ^"^ 1 ' 



iV-1 

X 

1=1 

(2.30) 



n(x(o) -A.)M(^-^)(Ai,--- ,A^_i)V^S"-^^(^')e-^^(-(o)). 



Prom (2.29), the resolvent in the 1-instanton background is given by 



1 / 1 
^tr 



_ 1 1 / 1 1 \ 

N"\z-(I)J/ ~ N z-X(o)'^\N z- Xi/ ■ ^^'^^^ 

\ »-l ' a;(o) 

As discussed below [sec (2.35)], this gives the resolvent in the 0-instanton sector at 
leading-order in the large- limit as 

1 / 1 \ \ /I / 1 \ \ (o-™st) 



=1 ' ^(0) 



(2.33) 

we see that the first term, 1/N ■ l/{z — ^(o)), corresponds to the 
isolated eigenvalue, which is distributed at x = X(o) ^ ^ (5-function, and the second 
term indicates that this eigenvalue causes a distortion of the distribution of the other 
eigenvalues, as shown by the factor n£i^(^(o) ~ Aj)^ in (2.30), both of which turn 
out to be at next-to-leading order. The physical meaning of is quite clear. 
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Because it is precisely the contribution from the instanton to the loop amplitude, it 
can also be computed by using the action (2.27) as surfaces with loops of q and q. 
Thus, it is found that AR^^~™^^^ describes the correction to the loop amplitude from 
open boundaries: 





It is not difficult to compute z\i?(^""^'^*) explicitly. First, we note that the re- 
solvent is also given by (2.7), even in the instanton background, because it can be 
derived generically using the Schwinger-Dyson equation. Thus, a different choice of 
the polynomial p{z) yields a description of a different instanton sector. Because the 
behavior of R{z) as z — > cxd fixes the coefficient of the first-order term in p(z), it is 
a constant term that distinguishes the instanton background. Denoting p{z) in the 
absence of the instanton by po{z) , the resolvent in the instanton background is thus 
given by 

1 



R{z) = - i^^V'{zY+pQ{z) + c + V'{z) 

= Ro{z) + - ^ = + • • • , (2.35) 

^,JV'{zY +Pq{z) 

where Rq{z) is the resolvent in the absence of the instanton, and c is a certain 
constant of 0{\/N) determined below. The above expansion corresponds to the 
right-hand side of (2.32) term by term. Therefore, AR^^~™^^\z) takes the form 

^^(i-inst)(^) = ^ I +.... (2.36) 

4y/y'(z)2 +Pq{z) 

In order to find the value of c in the case of the single-instanton background, we note 
that the eigenvalue density p{x) is related to R{z) as 

p[x) = --lraR{x + fO). (2.37) 
vr 

Using this, we impose the condition that the integration of Ap{x) corresponding to 
^^(i-mst)^^^ on the small interval around z = X(o) yields 1/A^, which implies that a 
single eigenvalue is located at z = X(o)- This condition amounts to requiring 
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where Cx^q^ is a small circle surrounding z = X(^oy This gives 

c=^W^(x(o)), (2.39) 

where 

Wo{z) = ^Vy^W+Mz)- (2.40) 
Substituting this into (2.36), we find 

In the double scaling limit, this quantity remains finite, and the first term gives 

8 w{C) 

as the instanton correction to the disk amplitude in the absence of the instanton 
w{C) given in (2.12). It is easy to check that we can obtain exactly the same result 
obtained from the case of the 0"^ potential. This implies that (2.42) is universal. In 
§4, we sec that (2.42) coincides with the results in the Liouvillc theory, or the loop 
equations. Therefore, we find that the instanton is identical to the D-instanton in 
c = noncritical string theory. 

§3. Chemical potential of the instanton 

We have seen that we can compute various quantities in the D-instanton back- 
ground by using the matrix model in the previous section. In this section, we show 
further that the matrix model also makes it possible to obtain a definite value of 
the chemical potential for the instanton (2.19), namely, the weight of the instanton 
itself. 

We start by computing the numerator in Eq. (2.19). For x < or X(2) < x, 
we can set det(x — cj)) = exp (tr ln(x — (j))) unambiguously. We hen obtain 



(det(x — (ff) = exp 



2 (trln(a;-<^)), + J ((2trln(x -</.))') + 



(3.1) 



■ 2 

where the subscript "c" denotes the connected Green function. Here, we have omitted 
tentatively the superscript "0-inst" . We note that 

(tr ln(x - </.)), = (tr ln(x - ct>))^,^^ + ^ (^) ' 
;(trln(a: - = <(tr ln(a: - <P)f)^^^^,^, + ^ (i^) ' ^^"^^ 
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where the subscripts "disk" and "cyhndcr" indicate the disk and cylinder amplitudes, 
respectively. Thus, Eq. (3.1) takes the form 



<det(x - 0)2) e-^^(^) 
= exp 



N (^2 ^ltrln(x-0)^^^^^-F(x)) +2 ((tr log(a; - c/.))^)^^^^^^^^ + O (^l) 

= exp (^-Nv}S\x) - V}i\x) + (^)) • (3-3) 

The 0(Ar)-contribution [-Nv}^\x)] consists of v}^\x) in Eq. (2.8) (as can be 
checked by differentiation), while the 0(A''°)-contribution V^^\x) is given by 

v}^\x) = -2 ((trln(x-<A))2) 



X — 



I cylinder 

.T(i) + ,T(2) 



-21n I 1 + y I +21n2. (3.4) 

^(X-X(i))(x-X(2)) 



In deriving Eq. (3.4), we have invoked the formula for the cylinder amplitude^^^ 

1 



tr ( tr 

x-^J \y- 



cylinder 



2{x-yy 



' xy- ^(x(i) +a;(2))(x + y) +X(i)X(2) 

^^(X - X(1))(X - X(2))(y - X(2)){y - X(i)) 



(3.5) 



Thus, the cylinder amplitude is given solely by the two endpoints x^X) and X(2) where 
X{x) < X(2), of the eigenvalue distribution, and it does not depend on any other details 
of the profile of the potential V{x). 

To evaluate the denominator of Eq. (2.19), it is also necessary to evaluate V^^^ (x) 
for X G X(2)]- However, Eq. (3.4) is invalid in this region. For instance, it 

diverges at a; = ^(i). This implies that we need to elaborate on the A/"-dependence 
of the cylinder amplitude. It is necessary to perform the calculation up to the 
magnitude of V^q , that is, the next leading order in If V^ff{x) in this region 
takes the form 

T/ W (x) = -lnN + [finite function of x] , (3.6) 

the chemical potential becomes finite, because the N dependence in Eq. (3.6) cancels 
the overall factor N in Eq. (2.19). The goal of the rest of this section is to demon- 
strate that this is indeed the case by evaluating the right-hand side of Eq. (2.19) 
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explicitly, using the orthogonal polynomial method. As we see below, {|det(x — (^)^) 
has a compact expression in terms of orthogonal polynomials. The task at hand 
then reduces to evaluating the explicit forms of orthogonal polynomials in the re- 
gions X £ [x(^i'j, a;(-2)], X < X(^i^ and x > X(^2) f'Cparatcly. It turns out that the case 
X G [2^(1), ^(2)] requires a more involved analysis, for which we employ a sort of 
WKB approximation explicitly (where h = jj) in solving the recursion relation for 
the orthogonal polynomials. By combining these ingredients, we finally obtain the 
chemical potential /i and show that is a "universal" quantity, independent of the 
details of the potential V{x). 

3.1. Orthogonal polynomials and (det(a3 — 0)^) 

It is known that a set of polynomials Pn{x) = + 0{x"'~^) (n G Z+ U {0}) that 
obeys the orthogonality condition 

J dx e-^^(^) Pn{x) Pm{x) = K Sum (3-7) 

can be expressed as the expectation value of det(a; — <^(„)) for a system of an n x n 
Hermitian matrix with the action ATtr 

Pn{x) = ^ / #W e-^*'-^('^w) det(x - , (3.8) 

where Z„ is the partition function of that system. Note the coefficient N multiplying 
the potential F((/>(„)) in the above expression. Equation (3.8) leads us to seek an 
expression of (det(a: — (f))^^ in terms of the orthogonal polynomials, for example, 
something like {P]si{x))^ . As a matter of fact, (det(x — 0)^) does have an explicit 
expression in terms of orthogonal polynomials and the coefficients of the recursion 
relation satisfied by them. We now demonstrate this point. 
Let us first define a new quantity Dn{x) as 



Z^,(a;)= det (M,,(x)), (3.9) 

0<i, ]<n—l 



where 



M,,(x)= / dx'e-^^(-') Y,^r^nx){x-Q)%Pk{x'). 

J irk yhihk 



(3.10) 

k=0 

Thus, Dn{x) is the determinant of the nxn matrix consisting of the first n x n entries 
of the matrix {Mij{x)}ij>o. In the above, Qij is defined through the recursion 
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relation for the orthogonal polynomials {-Pn}nez+u{o}) 

oo 

xPn{x) = ^ Qn,mPm{x) ■ (3.11) 
m=0 

In particular, for the n = N case we have 

L'iv(x) = (det(x - 4>f) . (3.12) 
It is easy to see that Dn{x) satisfies the recursion relation 

Dr,{x) = P^{X) + Qn,n-1 ^n-l(x) . (3.13) 

It is also straightforward to show that Qn,n+i = 1, while 

Qn,n—\ ~ fn ) 

r„ = --^(n>l), ro = 0. (3.14) 
Repeated use of the above recursion relation leads to 

Dn{x) = {PN{x)f + TN {PN-i{x)f + • • • + r^v • • • ri {Po{x)f . (3.15) 

This formula is the key equation in the evaluation of (det(x — <^)^), and it enables us 
to evaluate the chemical potential by using the asymptotic behavior of the orthogonal 
polynomials. The remaining task is to evaluate the orthogonal polynomials Pn{x) 
and the coefficients appearing in the recursion relation up to the next-to-leading 
order in the expansion. This is done in the following subsection. 

3.2. Asymptotic behavior of Pn 

The forms of the orthogonal polynomials in the large- A?^ limit have been deter- 
mined as follows:^^^'^^^ 

Pnix) = J- ef , ^ cos [N xa{x) -{N-n) Xb{x) + Xc{x)] • 

V ^ ^(x-a;(i))(x-a;(2)) 

(3.16) 

The values Xa,Xb ^'Hd Xc are evaluated in Refs.,^°) and.^-^^ Here, let us instead 
proceed to evaluate them directly from the recursion relation 

xPn{x) = Pn+l{x) + SnPn{x) + (3.17) 

where s„ = Qn,n in terms of Qn,m, in order to determine the behavior of the orthog- 
onal polynomials to next-to-leading order in O (-^). 
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At leading order in Pn{x) depends on A'^ as e^. Hence, it is convenient to 
consider the ratio 

M^) ^ - C?(A^°) , (3.18) 

and to expand the recursion relation (3.17) in With respect to kn{x), Eq. (3.17) 
takes the form 

T 

X = kn+l{x) +Sn + Trj-T ■ (3-19) 
rtn\X) 

We expand kn{x) in as 

Mx) = *'»'(.,;^)+ip*«(x.;^)+o(^). (3.20) 

When we take n and N to the infinity and introduce the continuum variable 
( = the leading and next-to-leading orders of the recursion relation (3.19) yield 

,^)y + {s{0-x)k^'^\x, e)+r(O = 0, 

, + (siO - x)) (x, + k^'\x, i\ disiO + d^k^^Hx, = , 

(3.21) 



X 



respectively. Here, the continuum limits of r and s are taken as expressed in Eq. 
(C.5) (see Appendix C). By solving these equations, we find 



(3.22) 



where the branch of the square root should be chosen so that kn{x) ~ a; as ^ oo 
by taking account of the original definition of kn{x). We also find 

The expression for A;W(x, 

can be further simplified as 

k^'H^, = —k^'^Hx, 09«ln«(^, 0, (3.24) 

where 



q{x, 0=2fcW(ar, ^) + s{0 - x 



(3.25) 
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Now, it is straightforward to write down the expression for Pn{x) as a product 
of the ki{x): 



Pnix) = JJ kj(x) = exp 



^lnfcj(a 



exp 



N 



Jo 



+ 



i(lnfc(°) (x, ^)-lnfcW(x, 0))+0(^ 



(3.26) 



where we have utihzed the following form of the Euler-Maclaurin summation formula 
to convert the summation into an integral 

n2 „!^2 



j=ni 



(3.27) 



Using the expression (3.24) for k^^\x) and the facts that r(0) = = s(0) and 
q{x,0) = x = A;(°)(x,0), we obtain an expression for the asymptotic behavior of the 
orthogonal polynomials as 



Pnix) = exp 



N I dClnfc(°)(a;, ^) + -lnA;(°) 
2 

1 ^ 
N 



n 



^' N 



+0 

'^(") (,x, ^) 



exp 



NT dClnfe(°)(a;, 

JO 



l + O 



(3.28) 



V 9 (a;, #) 

where k^^\x) and ^(a;, ^) are given explicitly by Eqs. (3.22) and (3.25), respectively. 

Thus, we have obtained an expression for Pn(a;), (3.28), for a given (real) value 
of X by solving this equation in terms of n. However, the above expression becomes 
complex when the argument of the square root in q{x, ^) is negative, i.e., when 

(x-.(^))=-4.Q<0. (3.29) 

This seems to be a contradiction, because Pn{x) should be a real polynomial. To 
resolve this difficulty, let us consider an actual plot of P„ [x) . There is a finite region in 
which it exhibits oscillating behavior with n nodes. This oscillatory region precisely 
corresponds to the region where the exponent of Pn{x) becomes imaginary in the 
expression (3.28), because Pn{x) contains q{x, ^) in its exponent. This suggests that 
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in this oscillatory region, Pn{x) in Eq. (3.28) should behave like a trigonometric 
function. 

Now, it is helpful to recall that wc have derived the expression (3.28) by starting 
from a set of linear difference equations. For this derivation, we can draw an analogy 
with the WKB method in usual quantum mechanics. The important condition here 
is the same as the continuity condition for the solution at the turning point, where 
the behavior of the wavefunction changes from decaying behavior in the classically 
forbidden region to the oscillating behavior in the classically allowed region. Thus, 
we should multiply the naive expression in Eq. (3.28) by an appropriate phase factor 
and take its real part. In other words, the phase factor emerging in the analytic 
continuation of Eq. (3.28) should be combined to make Pnix) a real trigonometric 
function with an appropriate phase. 

In Appendix D, we show that exactly the same argument as for the WKB ap- 
proximation applies here and, in particular, we should take an extra factor of 2 into 
account when Eq. (3.28) is continued into the oscillating region analytically. Noting 
also that \k^^\x,^)\ = ^/ri^, we obtain the following expression for Pn{x) in the 
oscillating region: 



Pn{x) 




exp 



NRe 



Jo 



sin {9{x)) 



(3.30) 



Here 6{x) represents a phase factor. As there are n 0{N) nodes located in a finite 
interval, 6{x) is expected to vary rapidly, while Eq. (3.30) explicitly shows that the 
amplitude of the oscillation changes very slowly. For our purposes, it suffices to know 
the averaged value taken over a small interval around x. Taking the average of the 
square of Pn{x) for an interval where 6{x) should oscillate at least once, we obtain 



(-fn(a;))averaged ~ ^ 



exp 



2iVRe / d^lnA;(°)(x, C) 





(3.31) 



which we use hereafter as the asymptotic form of the orthogonal polynomials when 
X lies in the oscillating region. 
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3.3. Evaluation of (^det(x — 0)^) 

Now we are ready to evaluate (|det(x — ^)^) using the asymptotic behavior of 
the orthogonal polynomials obtained in the previous section. There is an essential 
difference between the case in which x is positioned inside the cut, that is, < 
X < a;(2)5 ^-iid the case in which x is outside the cut. 

Let us first consider the case that x is located outside the cut. In (3.15), the 
ratio of neighboring terms, which we write as Cj{x), is 



Cj{x) = 



(Piv-j-i(x)) rN-j 



{PN-j{x)y 
r 



N-j 
N 



X 



N-j 
' N 

x) - 



+ 



{kN-j{x)y 
1 

N 



4r(0 



{s{0 -x) + ^(s(0-x)2-4r(0 



+ - 



c- N-j 
^ N 



(3.32) 



The quantity Cj{x) is of 0(iV°) and less than 1. Therefore, the largest contribution 
to (det(a; — 0)^) comes from the first term Pjv(x)^. Note also that the dependence of 
Cj{x) on j is modest [i.e. 0{N^)] compared to the number [0{N)] of terms. Hence, 
when the sum is replaced by a sum of terms with a constant ratio, say co{x), the 



correction will be of order 



AT 



<det(x - cl>f) = {PN{x)f J2 (^(^))' ( 1 + ^ 

3=0 

= {PN{x)f ^ 



1 - co{x) 



l + O 



(3.33) 



The factor 



l—co{x) 



turns out to be 



-(s(l) -x) + V(s(l)-x)2-4r(l)) 



1 - co(x) 2V(s(l)-x)2-4r(l) 
A;(o)(x, 1) 



q{x, 1) 



(3.34) 



Hence, we obtain 
(det(x - (t)f) - 



'fcW (x, 1) 
q{x, 1) 



exp 2N d^lnA;(°) (x, (l + O 



(3.35) 
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Next, we evaluate (^dct(a; — 0)^) for the case that x hes inside the cut, that is, in 
[x(x-), X(2)]- Note first, as briefly described in the previous subsection, the orthogonal 
polynomials Pn{x) have the following oscillatory behavior: (1) Pn{x) oscillates within 
a finite region in which it has n nodes. The amplitude of each oscillation itself 
varies slowly, so that we can draw a curve enveloping these oscillations. Outside 
this oscillatory region, Pn{x) behaves as (2) This oscillatory region grows as n 
becomes large. For sufficiently large n, this region becomes the region in which k^^^ 
becomes imaginary, 

(.-»(^))^-4.(^)<0. (3.36) 

Thus, the main contribution to Dn comes from the terms contains Pm{x) with 
m> HQ, where the integer no for a given x is defined by 

.-M=f^))^-4.(^)^0. (3.37) 

We can repeat the argument that leads to (3.33) for this case. We can thereby 
show that the contribution from Pm{x) for m < no{x) is suppressed exponentially, 
because the ratio of the neighboring terms is less than 1. Hence, (det(x — (j))'^) can 
be approximated as 

N N-l-i 

{det{x - cl>f) = J2 (^^(^))' n ^^-J - 

i=no{x) j=0 

Prom the expression for in (3.31) and an additional calculation presented in 

Appendix E, we obtain 



<det(x - ^)'>,,,,,g,d = 2^^ Pi^) 



2N f ° d^\nk^^\^) + N [ d^lnr(0 
Jo ho{x) 



l + (3.39) 

where ^0(2^) = and p{x) is the eigenvalue density (2.37). 

We summarize the above result in the form of the effective potential V^fj from 
the relation e-^(^«ff (^)-^(^)) = (det(a; - (t>f). 

Outside the cut: V^^{x) = V{x) -2 [ d^lnk^^^O 

Jo 

-l{-log[(x-5(l))2-4r(l)] (3.40) 
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+2 log (x - s(l) + ^/{x - s(l))2 -4r(l)) - 21og2|. 
Inside the cut: Ves{x) = V{x) - 2 d^lnk^^' {^) - d^lnr(0 



-ilog[2Ar7rp(x)Vr(l) 



(3.41) 



Also, if we use V^q in Eq. (2.8), the above expressions can be cast into much simpler 
forms. 

'A;(0)(x,1) 



Outside the cut: Ves{x) = v}^^ ~ ^ { ~^ 



q{x,l) 



Inside the cut: Veif (a;) = v}q\x) - log 2N^^p{x)^/r(T) 



(3.42) 
(3.43) 



3.4. Universality of the chemical potential of the instanton 

Now we are ready to take the ratio of the partition functions existing inside 
and outside the cut. They are given by the integral of Ves above obtained for the 
respective regions. We evaluate this ratio, the chemical potential of the instanton, 
near the critical point. 

We can evaluate the contribution from outside the cut using the saddle point 
method. To do this, let us first evaluate the leading contribution to VeS, that is the 
first two terms on the right-hand side of (3.40), 



V{x)-2 /'deinfe(0)(O. 
Jo 

ppendix C, r 



Near the critical point, as explained in Appendix C, r and s behave as 

1 



r = Tr avr 



(3.44) 
(3.45) 



where Ag represents the deviation from the critical point and a is a certain constant. 
If we rescale x and ^ to emphasize the region near the critical point as 

x = .s, + 2V?;+ V^C, (3.46) 
^ = 1 - (3.47) 

Vf^s can be expanded in terms of Ag as 



Kff ~ V{x) - d7]Ag I 2 ln( V^) + Ag^ ■ 



1 2a/ Cv^ + rcOiVrj + 1 



+ OiAgr 



(3.48) 
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We arc interested in the universal part of the derivative of Vqs rather than Vcs 
itself. In other words, we arc interested in the term proportional to Ag^ , which can 
be obtained in the following familiar form after integrating over ij: 

y^^^^^_S i(SziM<±^, (3.49) 

Thus, at leading order, the integral / dxe''^"^ can be evaluated at the saddle point, 
C = sv^*^- next-to- leading order, that is, for the contribution from the 

cylinder part, we can use the same saddle point as that at leading order. The 
standard saddle point calculation yields 

Jx<x^i-),x^2)<x l^ ai Ag8-\/N 

The contribution from inside the cut is simpler to evaluate. Noting that for this 
case, V^g{x) takes a constant value, say V^g{x(^2))i we can replace rjv with Tc, 
it is straightforward to obtain 

Jx^i)<X<X^2) Jx^i)<X<X^2) 

(3.51) 

In passing, we note that we can perform another calculation that leads rather directly 
to (3.51). This calculation, which we present in Appendix E, does not involve the 
explicit form of V^fj(x). 

Taking the ratio of (3.50) and (3.51), we finally obtain the chemical potential of 
the instanton (2.19) as 



3424 _jv(yW(a.(o))-yW(x(2)))_ ^3_^2) 



To make a connection with the standard analysis, we need to rewrite the chemical 
potential (3.52) in terms of a quantity that appears in the string equation. The 
solution of the string equation is the second derivative of the free energy, and it 
is conventionally normalized so that F ~ —\/t. With this normalization, the free 
energy is given by 

-At§. (3.53) 
On the other hand, the free energy of the matrix model is obtained as 

i^matrix = f\l - ^og{r{0)dC (3.54) 
Jo 
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Near the critical point, r behaves as (3.44). Expanding i^matrix in terms of Ag and 
taking the first non-integer exponent yields the universal part of the free energy as 

.^J^. (3.55) 
15 

By comparing (3.53) and (3.55), we can connect Ag with t in the string equation as 

2 

Ag = ^. (3.56) 

Using this relation and invoking Eq. (2.10), we reach the universal expression for 
the chemical potential 

^(l-inst) . _ _ 

1^ = = 3 re"~* . (3.57 

If V{x) is an even function, the above argument needs two modifications, due 
to the accidental symmetry of the potential with respect to the exchange of x and 
—X. First, s vanishes identically. However, (3.48), (3.49), (3.50) and (3.52) hold if 
we replace a with ^ in these equations. Then (3.52) becomes 



1 1 

3422 

:e 



^{v}sHm)-y}sM)_ (3.58) 



Secondly, there emerge two critical points corresponding to two maximum points of 
the potential that have exactly the same height. Each critical point contributes the 
same amount to the free energy, and therefore the value of the free energy is twice 
that in generic potential cases. Taking this effect into account, (3.53) should be 

8 5 , , 

--t-2. (3.59) 

On the other hand, because s docs not contribute to the free energy, (3.55) does not 
change. Thus the relation between Ag and t becomes 

4 

^9 = 4^- (3-60) 

Substituting this relation into (3.58), we reach the same result. Hence, it is concluded 
that (3.57) is universal. As a concrete example, we present a computation of /x in 
the cases of both the 0^ and 0^ potentials in Appendix F. 
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§4. D-instanton effect in loop equations 

The Schwinger-Dyson equations for the correlation functions of the "loop vari- 
ables" tr^" are known as the loop equations, and it is believed that they give a 
complete description of the system. Their continuum limit is easily taken,^^) and 
the result can be interpreted as a sort of closed string field theory. The loop 
variables are also useful when we compare the matrix model with the Liouville field 
theory. 

In this section, we consider the question of what properties of the instantons 
observed in the previous sections can be captured by the continuum loop equations. 
After some preparation in §4.1, we examine the classical limit (^st 0) of the 
loop equations in §4.2, and we see that in the single-instanton vacuum, the classical 
solution for the loop amplitude has an ambiguity in the form of an arbitrary constant 
factor. In §4.3, we show that this constant factor can be determined if we consider 
the loop equations to all orders. In §4.4, we first show that the single-instanton 
loop amplitude obtained in this manner indeed is identical to that obtained from 
the matrix model calculation. We next compute the D-instanton effect on the loop 
amplitude in the Liouville field theory, and we find that it too reproduces the result 
of the matrix model. This coincidence confirms that the instanton in the matrix 
model is identical to the D-instanton in the Liouville field theory. 

To this point, we have seen that the loop equations correctly describe the loop 
amplitude in each vacuum with a fixed number of instantons. However, the con- 
tinuum loop equations cannot determine the chemical potential of the instanton. 
In §4.5, we see that the continuum loop equations give a divergent expression and 
require regularization when we attempt to evaluate the chemical potential. 
4.1. Loop equations 

In this subsection, wc make a few remarks on the loop equations that will be 
important in the succeeding analysis. If we start with the matrix model with the 
action 



S = NtTV{(l)), 
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the loop equations can be derived from the Schwinger-Dyson equations 

= y # ^ [tr Wn,--- Wn^e-^''^^t>'>] . (4.2) 

Here, W„ is the loop variable 

T^„ = ^tr</,-, (4.3) 
and (j) = (l>""t"' with the U{N) generator normalized so that ivitH^) = S"-^. Using 

tr(t"^) tr(r5) = tr(^5) 

tr (i"A fB) = tiA tvB , (4.4) 
we can show that Eq. (4.2) takes the form 

no-l p 

{WjWno-l-jWn, ■■■Wn^)-Y^ Cm {Wn,+m-lWn, ■ ■ ■ 
j=0 m=l 

+ ^2 E (^-1 • • • ^n._i " " " ^n.) = . (4.5) 

s=l 

We would like to find the form of Eq. (4.5) in the double scaling limit, expressed by 

1 5 ^ ^ 

— = a2 , (4.6) 
for ^^-theory in which V{x) = ^ — | x^, and by 

™ _ ™ pa7C 

1 5 

j^=a^ i^9s, (4.7) 

for (?i)^-theory in which V{x) = ^ — f a^"^- In Eqs. (4.6) and (4.7), and g^, represent 
the critical values of x and g, whose values, as well as those of 7 and k, are given 
explicitly in Appendices A and B for 0^-theory and (^^-theory, respectively. If we 
take the continuum loop length I as I = ^na, Wn approaches the continuum loop 
operator 'w{l) in the double scaling limit as 

Wn^Jx^+'jZlwil), (4.8) 
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where the constant Z^/ is given by 

Zw = ^r'{g^x^f for (/.^-theory (4.9) 

and 2 

Z^^=(^ 7^ for 0^-theory . (4.10) 

Using this, we can show that the loop equations (4.5) have the fohowing continuum 
form in the double scaling limit 

I f dl' {w{l') W(l - I') W(li) ■ ■ ■ W(lr)) + Psil) {wih) ■ ■ ■ W{lr)) 

Jo 

r 

+Sst I XI ^« ^^'^^l) ■ ■ ■ '^Q-s-l)wils + l)w{ls+i) ■ ■ ■ Wilr)) = , 
s=l 

p,il)^3S"il)-ltdil). (4.11) 

We note that ^st appearing in Eq. (4.11) differs from gg in the matrix model by 
the factor 

1 K 

gst = -= gs. (4.12) 

■s/Zw ix^ 

Using Zw in Eq. (4.10) and the values for /3, 7 and n given in Appendices A and B, 
we find that th 
and <j)^-ih.eoi:y, 



we find that the conversion factor ^1— takes the same value in both 0^-theory 



We also remark that we employ such a convention that the same renormalized cos- 
mological constant t and renormalized boundary cosmological constant C, can be used 
in the loop equations and the matrix model. Therefore, we can compare the quan- 
tities in the loop equations with those in the matrix model by taking account of the 
conversion factor (4.13) for the string coupling constants. 

To treat the coupled equations (4.11), we introduce canonical pairs ^(Z) and 
ip^l) of closed string fields that satisfy the relations 

'm.i'Hi')] =s{i-i'), 

[V(Z), VOO] =0 = [v^(0,V'^(O] , (4.14) 

and the "vacuum" state |0) that is annihilated by all '^(Z) i.e., for which we have 
^(/) |0) = ((0| VUO = 0)- We then define the state by 

l'^) = E -| / V'H^l) • • • / dlr i^^lr) |0) {wiW) ■ ■ ■ wilr)) 

r=0 
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exp 



°° -1 POO /"OO 

V- / dh^\h)--- dlri^^{lr)w^'Hh,--- ,lr) 

,r=l •'^ •'^ 



|0) , (4.15) 



where w^"^^ (hr ' ' > represents the connected part of the Green function {w(li) ■ ■ ■ w{l, 

(4.17) 

Using \^), the original Green functions can be expressed as 

{wih) ■ ■ ■ Wilr)) = (0| i^ih) ■ ■ ■ ^{Ir) W) , (4.18) 

and the continuum loop equations (4.11) can be written in the compact form 

D{1)\^)=0, (4.19) 

where 

Dil) = ITil) + p,il) , 

rl poo 

T{1) = / dl' I') + g'i / dl' I' {I') + /') . (4.20) 

Jo Jo 

For subsequent considerations, we remark here that T{1) is isomorphic to the 

energy-momentum tensor in momentum space if I is regarded as a momentum. To 

see this, we construct a new variable ip{l) from 1^(1) and V'UO ^ 

r V2 



m = { 



9st 

ffst 

I V2 



■m 



(/>0) 



(4.21) 



(-0V^(-0 (Ko) 



This variable satisfies the commutation relation 

[(^(0,(^(0] = m+i'). 

T(l) can then be rewritten as 

/OO 1 
dl' -■.^{l')ip{l-l'): 
-OO ^ 

«;W(«i,---,/r)=0((ff.tf('"-'^) , 



We recall that 



(4.22) 

(4.23) 
(4.16) 



in our convention. 
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If we further introduce their "coordinate" representation by 

^ J —oo * 

/oo 

dl e'^" T{1) , (4.24) 
-oo 

r(cr) can be expressed as 

TicT)=gll:{daip{a)f:, (4.25) 

and the two-point Green function of da-<p{cr) is given by 

(0| 5,^(a)5,,^(a') |0) = ^ , . (4.26) 

(a - a') 

Hence, T'(cr) is the energy-momentum tensor of a free massless scalar field ^(cr) up 
to an overall factor of g^^. 

4.2. Classical approximation of the loop equations 

This subsection deals with the classical limit, gst —>■ 0, of the loop equations. In 
the perturbative expansion, the amplitude w^^\l) of one external loop begins with 
the disk contribution w^\l), which is 0{{gst)^) in our convention (4.16). Because 

u;(2)(/^^ I2) = 0(5^4) we have 

(^(/i)^(/2)) = wi^\k)wi^\l2)+0{gi) . (4.27) 

Thus, Eq. (4.11) with r = at C((5st)°) becomes 

d^{w{C)f + 0{{g,^f) = 3e-lt, (4.28) 

where w{C) is the Laplace transform of w^^\l). This equation determines {w{Q)'^ 
up to a constant zu: 

(^5o(C))' = (C - ^) (C + Vt) . (4.29) 

In the classical limit, the loop equations do not place a restriction on w. 

As we have seen in §2, in the matrix model, w is determined by the number 
of eigenvalues positioned on the top of the potential. By contrast, in the classical 
approximation of the loop equations, w is ambiguous. This is in some senses obvious 
from the beginning. It is clear that w cannot be described in terms of classical 
solutions of closed string field, because it is not of order g^^ but rather ^st- 
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4.3. Loop equations to all orders and the nonperturbative effect 

Beginning in this subsection, we investigate the possibiUty of determining the 
constant w and the chemical potential using the loop equations to all orders, (4. 19). 

First, we examine the vacuum expectation value (VEV) of an operator A in 
the presence of a fixed number of instantons. This means that we pick up, say, a 
single-instanton sector out of an infinite number of sectors and consider the VEV of 
A in this sector. In the matrix model, such a quantity can be written as 



{A) 



(l-inst) 



'N-1 



(0-inst) 



X N 



I 

J Xr- 



dx {A{^\ x) det(x - e-^^(-) 



"=(0) 



cx iV / dx {A{4), x) exp [2tr ln(x - 0')] ) 
= N [ dx e-^^(-) 



/(0-inst) Nv(x) 



X ( A{4>', x) exp 



-2 / dl 

10 



e 



d 



tr e 



14>' 



— tr 1 



/ (0-inst) 

(4.30) 



where Z 



I (0-inst) 
JV-1 



and (, ) 



/ (0-inst) 



are defined in Eq. (2.16). 



Because we instead analyze the loop equations here, we need to construct a state 
describing the single-instanton vacuum by closely considering the expression (4.30). 
Equation (4.30) expresses the VEV of an operator A in the single-instanton vacuum 
in terms of the VEV in the null instanton vacuum of the operator obtained by making 



the loop of length I (tr e}'^' ) condensed with a weight 



and integrating with 



respect to x around a;(o)- In the continuum limit, the null instanton vacuum state 
l'^)o-inst is given by 



1^^)0-1 



mst 



V - / dh ■■■ dlr ^\lr) W^'^ ih,--- ,lr 

,r=l "'0 



exp } ^— I (Ui1p'{li)--- I (Ur1p'{lr) W^'' [h,--- Jr) , |0), 

(4.31) 

where w^'^\li,--- , ^r-)|pgrt w^'^\li,--- , Z^) in the null instanton vacuum. The 
single-instanton vacuum state can be obtained by applying the operator V(C) to 
l'^)o-inst' which condenses the loops with an appropriate weight, and then integrat- 
ing V(C) 1^^) 0-inst '^^^^ respect to C around ^ , corresponding to X(o)- The operator 
V(C) is obtained by replacing the operator -^tre'*^' in the exponent of Eq. (4.30) 
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with 



V(C) = C'inst exp 



-Q 



' 2 



(4.32) 



Here, we have introduced the parameter Q, to take account of a possible renormal- 
ization of the order iV ~ and the overall normalization Cmst- Using V{Q in 
Eq. (4.32), we obtain the single- instanton vacuum in the form 



inst 



2 



= a 



inst 



^ 2 



exp 



inst 



-Q 



2 

y e 



inst 



(4.33) 



Our investigation to this point has not considered the all-order loop equations 
(4.19). We explain in the following subsections that Q is identical to the constant 
ambiguity w in the loop amplitude and that Cinst is related to the chemical potential. 
In the rest of this subsection, we would like to find whether Q and Cinst can be 
determined by analyzing the all-order loop equations. 

As we have seen above, the operator V{Q in Eq. (4.32) has been proposed 
to introduce an additional instanton into the system if ^ is integrated around 



Prom the viewpoint of the matrix model, 



^ 2 



dCV(C) plays the role of adding 



an eigenvalue at the top of the potential. Hence, / dC,V{C,) adds an eigenvalue 

J —oo 

without specifying its position. We recall that the full vacuum consists of the 
vacua of various numbers of instantons: 



1"^) = l'^)o-inst + l*')l-inst + l*')2-inst + 



Thus, we should have 



(iV+l) 



/oo 
-oo 



(AT) 



(4.34) 



(4.35) 



where \^)'^^^ is the vacuum of the system consisting of N eigenvalues. In the large- 
N Jimit, this equation requires that if \^) is a solution of the loop equations (4.19), 
V(C) 1^) must also be a solution of these equations. This will be the case if 



^(C),v(C) ~WC) 



(4.36) 



holds, where D{Q is the Laplace transform of D(l) . As remarked in §4.1, the 
operator part of -D(C) is essentially the energy-momentum tensor T{Q. Thus, the 
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above relation implies that V(C) should be a primary field of conformal dimension 
equal to 1. The works in Ref.^^) carry out a search for such an operator V(C) and 
find an expression in terms of local operators. We find that the primary operator 
with conformal dimension 1 that resembles Eq. (4.33) is 



: exp ( —iV2 (p(a)] : = : exp — / dl ( — - > 



-ila 



^t(Z) 



(4.37) 



To bring this into the form of Eq. (4.33), we apply the analytic continuation a — ^ i(, 
which yields 

2 



V(C) = Cinst : exp 



/ 

Jo 



dl 



gstl 



(4.38) 



We thus find that the single-instanton vacuum |'?^)i_iiist given by 

|'^)l-inst=Anstl'^) 



/0-inst ' 



D: 



inst 



Cinst / ^dC -.exp / 
Vc~^ Uo 



dli ,e-^^V(0+5ste^^V^(0 



9stl 



. (4.39) 



The expression (4.39) shows that the all-order loop equations determine Q {Q 

/oo 
dC V(C) with V(C) in Eq. (4.38) is 
-oo 

always a solution of the loop equations for any value of Cnst- Thus, the normalization 
constant Cnst cannot be determined as long as we use only closed string fields. 
4.4. Loop amplitude in closed string field theory 

Now we know how to express the instanton contributions in the loop equation 
formalism, at least up to an overall normalization. In this subsection, we exam- 
ine how they appear in the limit ^st 0- Then, we compare the results of the loop 
equation with those of the matrix model and the Liouville theory approaches. In par- 
ticular, as we discussed in §4.2, they should correspond to the solution of Eq. (4.29) 
with some w. We calculate the value of m using the loop equation and compare it 
with that obtained in the other approaches. 

In order to study the limit ^st 0, it is convenient to express the state (4.39) 
using Eq. (4.15) as follows: 



l'^)i-i 



inst 



c 



inst 



^ 2 



X exp 



d(^ exp 



oo ^ 
r=l " 



dlgste'^'i^^l) 



dh (i^\li 



gsth 
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-CI, 



pert 



|0). 



(4.40) 



Because w^'''\li, ■ ■ ■ ,lr) = O ((gst)^'-'^^^'') , we can expand Eq.(4.40) in terms of g^t- 
Using this expression, let us calculate {^\^)i-i^st ^'^^ (^1 ^(0 l'^)i-inst' which can be 
compared with the quantities calculated in the matrix model formulation. 

The quantity (0|^')i_jjjs|- should be the continuum version of j^^^^^^^ (ix 6^^°**^^^ 
in the matrix model, and it can be expressed as f ^ d( e~^''^^'^^ in terms of the 
continuum effective potential Ves{C)- Thus, using Eq. (4.40), we find that V"eff(C) 
should be expanded as 



dl w^^\l) 





1 

+ 2 
+ •• 



pert 

2 



-Cl 



CO 

dh , 

S'st'l 



oo o 

ysth 



pert 



(4.41) 



The leading-order contribution to w'^^^l) is from the disk amplitude w^\l). Thus, 
the above formula to leading order coincides with the matrix model result in the 
limit AT — > oo. The leading-order contribution to the second term on the right- 

(2) 

hand side is from the cylinder amplitude, Wq (Zi, Z2). Higher-order terms involve 
the contributions from worldsheets with more boundaries. Therefore, the expansion 
above corresponds to the expansion Eq. (2.28) in the matrix model. 

The quantity {0\'tp{l) |!P')i.inst should correspond to the loop amplitude in the 
instanton background w{l)\^^^. It is calculated as 



= Qiist 



Jc-^ 

Cl 



gste 



°° 1 / \ / 



-Clr 



9stl 



r-1 



X W^''\li,--- Jr-1, I) 



pert 



(4.42) 
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The first few terms in tfie gst expansion are given as 



pert 



(4.43) 



In tlie limit g^t 0, we can use the saddle point approximation to evaluate the 
integration over C,. The saddle point Q = \/i/2 can be identified with X(o) i'^ the 
matrix model. Thus, by Laplace transforming the above expression, it is easy to see 
that the three terms in it coincide with the matrix model results Eqs. (2.32) and 
(2.33). 

From the analysis to this point, it is obvious that what we have been doing 
here is exactly the continuum version of what we did in §2. To be more concrete, 
Eq. (4.35) should be considered as the continuum version of the manipulation used 
to isolate the contribution of one eigenvalue by introducing qi and qi. The value C, 
in Eq. (4.35) corresponds to the eigenvalue x in the continuum limit. 

The cylinder amplitude Wq {li, I2) can be obtained by solving Eq. (4.11). It is 
given as 

4'^ (Ci, C2) = -gl In ijCi + ^t + \lc2 + Vt^ (4.44) 

in the Laplace transformed form. Thus '5;|inst(C) = dle~'^^w\i^t{l) is obtained as 

= ^'^^0 + l'-^ + Oiig.,)^). (4.45) 
By comparing this with the expansion of w{Q in zc given by Eq. (4.29), 

^5(0 = MO + I + 0((5st)') , (4.46) 

we find that 

^lloop equation = 5st ^6 . (4.47) 

Taking the factor in Eq. (4.13) into account, this value of zu is found to coincide 
with the matrix model result (2.42): 

^Latrix = 5.^i^- (4.48) 
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This value of w is consistent with the result of the Liouville theory. In the Liou- 
villc theory,^^) the instanton wc have been studying corresponds to the D-instanton. 
The amplitudes in the presence of such a D-brane can be calculated by using the 
open string theory. In particular, ttj(Z)|inst can be evaluated as an expansion with 
respect to the string coupling constant. The first correction to the disk amplitude 
is given by the cylinder amplitude with one boundary on the D-branc. Thus w can 
be obtained by calculating such a cylinder amplitude. This is exactly what we have 
done above. 

In the Liouville theory,^'') the D-instanton corresponds to the ZZ-branc bound- 
ary state with (m, n) = (1, 1). The loop w{l) that we have been dis- 
cussing corresponds to the FZZT boundary state |-Bs)fzzt- normalization of 
\^{m,n))iz ^'^^ I^s)fzzt fixed through a modular bootstrap or, in other 
words, the open-closed duality.^) Therefore, we need to calculate the cylinder am- 
plitude with one boundary on the ZZ-branc and the other on the FZZT-brane. This 
can be calculated in a manner similar to that employed in Ref.,^"*^^ but here we pro- 
ceed differently, as we now demonstrate. It is known that the properly normalized 
boundary states jj^^^z I-^s)fzzt satisfy the following relation^^^'^^^'^^) : 



l^{m,n)>zz = ^.=i(f +n6))p22T " ^^=<f -«^) /pZZT " ^^"^^^ 

We have h = -y^ for the c = case, and s is related to the boundary cosmological 
constant jiB and the bulk cosmological constant jj, in the Liouville theory as 

cosh(7r6s) = (sin(7r62))^ . (4.50) 

The relation (4.49) can be used to calculate the cylinder amplitude in question from 
Wq\Ii, I2) obtained using the loop equation. Thus, we can check if our results are 
consistent with those obtained from the Liouville approach. 

The value of jiB-, or equivalently s, is obviously related to the variable C, we 
have been using. Because h = s = z (| it 6) gives cosh(7r6s) = |. On the 

other hand, the Laplace transformation of the disk amplitude wq is proportional to 
cosh(7rs/6). Thus we see that wo vanishes if s = i (^ ± 6). Therefore s = z (^ it 6) 
should correspond to C = • A careful analysis shows that s = i[^ + b) corresponds 
to C = on the second Riemann sheet, while s = z (^ — 6) corresponds to that on 
the first Riemann sheet. 
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Now we calculate the cylinder amplitude. The FZZT-brane boundary state can 
be naturally identified with the boundary formed by — trln(2; — cp) in the double- 
scaled matrix model including the normalization. Therefore, the 0{gst) term in 
w{Q should be 



dl e-^' 











T 



9st 



on the second Riemann sheet 



C'=^ on the first Riemann sheet 



(4.51) 



Substituting the cylinder amplitude given by Eq. (4.44) into Eq. (4.51), we obtain 

1 



9^ 

1 



f' = ^ on the second Riemann sheet 



5st 



1 



on the first Riemann sheet 



2 Vc+Tt Vv/c+7t-^/iVi v/c+7t+^/i^, 

1 ffstVsti 

2 w^{Q 



(4.52) 



This shows that vd = gstV^t^, as anticipated. 

The relation in Eq. (4.49) can be seen more directly from the matrix model 
point of view. Specifically, as we have seen in Eq. (4.30), adding a single boundary 
corresponding to the ZZ-brane amounts to inserting 2Retrln(2; — (f)) evaluated at 
z = X(o) from the matrix model point of view: 



l%i)) 



zz 



2Retrln(z 



(4.53) 



On the other hand, as we noted above, the FZZT brane corresponds to the integral 
of the resolvent from the first sheet to the second sheet: 



B 



FZZT 



B 



FZZT 



C{0) 



on the second sheet 



dz tr- 



t=(()) 



on the first sheet 



2Retrln(z 



1^ = ^(0) 



(4.54) 



Because the relation (4.49) was obtained by using the open-closed duality, Eqs. (4.53) 
and (4.54) suggest that the information concerning the open-closed duality is some- 
how incorporated in the matrix model. 

4.5. Ambiguity in the normalization of single-instanton vacuum state 

As we have seen in §4.3, we cannot determine the constant Cinst using even the 
all-order loop equation. For this reason, we cannot determine the chemical potential 
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fi using this approach. However, as we have shown in §3, /i ean be determined as a 
universal quantity in the matrix model approach. Therefore, if in the loop equation 
approach we can mimic as closely as possible the procedure applied in the case of 
matrix model, we should be able to calculate /i. As we argued in the previous 
subsection, Eq. (4.35) can be considered the continuum version of the manipulation 
used to pick up the contribution of one eigenvalue in the matrix model. Hence, we 
can use this fact as a guide to obtain the continuum version of the calculation of the 
chemical potential. Although is cx) in the double scaling limit, here we consider 
A^ to be large but finite and take the limit A^ — > oo later. Dividing the integration 
region into two parts as 

=( r dCV(C) + / ^'dCVio] (4.55) 

\J-y/t J-oo I 

in Eq. (4.35) (here we have replaced A^ by A^— 1 for consistency with the matrix model 
formulation) , we can regard the second term on the right-hand side as corresponding 
to the zero-instanton sector and the first term to the single-instanton sector. In this 
formulation, \i can be derived from the ratio of the first term to the second term, and 
it does not depend on the overall normalization of the operator V. Let us examine 
whether we can calculate /x in this formulation. 

In order to calculate /U, we should evaluate (0|ll')^'^\ which can be expressed 
as (iC 6~^<=ff(''). It is very easy to obtain identities similar to Eqs. (4.40) and 
(4.41) in this case. Because C, > —Vi corresponds to the single-instanton sector and 
—\/i>C does to the zero-instanton sector, we calculate V^fr(C) separately for these 
two cases. 

For C > -Vi, we can use Eq. (4.41) without any change. In the limit Qst 0, the 
leading contribution is from the disk amplitude, and it reproduces the matrix model 
result in the N ^ oo limit, up to an additive constant. The cylinder contribution to 
—Ves{() is obtained as 

-21n (^2^C + Vi^ . (4.56) 

We expect that the cylinder contribution, which corresponds to the one-loop ampli- 
tude for an open string, should contribute to the chemical potential fi. 

For —\/t> C, in evaluating (0|<P')'^^\ we should fix the value of the loop ampli- 
tudes on the cut in the complex C plane. For our purposes, here we should mimic 
the manner in which we dealt with such an ambiguity in the matrix model. Doing 
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f oo o /"OO 



this, we obtain 

/•oo 

-^efr(C) = Re^/ dlw^^\l) 
Jo 

-I /-oo 9 /•oo 2 

+- Re^i / d/i — - e-^i^i Rec, / d/2 —r e"^^'' (^i, ^2) 
^ 7o 5st'i 7o 5st'2 Ci=C2=C 

+ • • • . (4.57) 

Here, Re^ denotes the operation of taking the real part as a function of (. What 
is relevant for us is the difference V^fr(^) — Vefl{—^/i) and that the leading-order 
contribution precisely reproduces the matrix model result. The second term can be 
evaluated as 

- ^ In (^\/Ci + Vi + \JC2 + Vi^ + In (^yj (i + Vi - \fc2 + Vi^ 

+ln (^-^Ci + Vi+ v^C2 + V^) + In (^-\/Ci + Vi- \f^+Vt^ 

= -ln(Ci-C2)lci=C2=0 (4-58) 

which is divergent. The cylinder contribution is, in a sense, continuous at C = 
— because Eq. (4.56) is also divergent when C, — —\/t. In any case, because of 
this divergence, we cannot reproduce the value of the chemical potential ji in this 
approach. 

If we proceed ignoring this divergence, we obtain the chemical potential as 



t~^-l X N X oo"\ (4.59) 



where the factor t'~2 comes from the cylinder contribution in Eq. (4.56), from 
the Gaussian integration around the saddle point, and N from the combinatorial 
factor. Therefore, the loop equation approach reproduces the essential part of ^, i.e. 
t~ 8 , but it fails to reproduce the precise numerical factor. What we have seen in this 
subsection is the continuum version of what we observed in the first part of §3. In 
the matrix model, the divergence is somehow regularized in conjunction with N, and 
we obtain a finite value for fx. In order to calculate fj, in the continuum approach, 
we might need some renormalization procedure. 

Below, we summarize the results obtained for the nonperturbative properties of 
a noncritical string via the loop equations, or equivalently, the closed string field 
theory: 
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(1) In the classical approximation (gg 0) of the loop equations, in a background 
of D-instantons, we can obtain information regarding the D-instantons through 
the parameter Q in the instanton creation operator (4.33), which is related to 
the constant zu in the loop amplitude (4.29). However, the value of Q cannot be 
determined within this approximation. 

(2) The treatment to all orders of the loop equations determines Q. 

(3) The chemical potential of the instanton cannot be determined using the loop 
equation approach. 

It may be useful to compare these three results with the results concerning D-branes 
obtained in critical string theory: 

(1) The boundary state for the D-brane is constructed in the closed string theory. 
The normalization of the boundary state cannot be fixed if only closed strings are 
considered. 

(2) The normalization of the boundary state of the D-brane is determined through 
comparison with the one-loop amplitude of the open string. 

(3) In order to determine the probability of the creation or annihilation of D-branes, 
which is the same as the chemical potential /x, we need knowledge concerning the 
nonperturbative dynamics of an open string.^'')' "^^^'^^^ 

Although we do not know the critical string version of what we have done for non- 
critical string theory using the loop equation approach, we can see that there is a 
similarity between the situation for critical string theory and that for noncritical 
string theory. 

§5. Conclusion 

One possible non-perturbative formulation of string theory is that which uses 
closed string field theory or loop equations. The situation would be simple if string 
theory could be formulated non-perturbatively using only closed strings. However, 
what is implied by the facts we have just determined is that a formulation based on 
only closed strings may not be able to incorporate all of the non-perturbative effects. 
For example, to calculate the chemical potential for the instanton, we need to know 
the cut-off dependence of the cylinder contribution. As we have seen, this dependence 
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is well-defined in the matrix model, but not in the loop equations formulated in 
continuum variables. 

Using the loop equations in the classical limit, that is, the large- limit, one can 
determine the loop amplitudes in an instanton background, up to a constant factor, 
which is denoted by Q in §4. Treating the loop equations to all orders, one can 
further determine this factor. In terms of the D-brane of critical string theory, these 
two procedures, respectively, correspond to determining the boundary state up to 
the normalization from the boundary condition and determining the normalization 
of the boundary state using the one-loop calculation of the open string in the dual 
channel. Thus, the calculation to all orders of the loop equations is powerful enough 
to determine the factor that corresponds to the normalization of the boundary state, 
without the need for the calculation in the dual channel. Despite this fact, however, 
it cannot determine the chemical potential for an instanton. This corresponds to, 
in critical string theory, determining the weight of the annihilation and creation of 
D-branes from the dynamics of the tachyon. We have shown that the matrix model 
does determine the chemical potential universally. Thus, at least for noncritical string 
theory, the closed string does not describe the nonperturbative effects completely and 
it seems that the description offered by the matrix model is more fundamental. If this 
is also the case for critical string theory, it strongly suggests that the non-perturbative 
formulation of string theory must contain degrees of freedom corresponding to open 
strings or matrices. 

Because the chemical potential of the D-instanton is a universal quantity, it 
is conceivable that there exists some continuum approach to calculate it. Open 
string field theories for noncritical strings have been constructed in Refs.^^^-.^^^ It 
is an intriguing problem to calculate the value of the chemical potential using such 
theories. 

The calculation of the chemical potential employed in this paper is also applicable 
to other matrix models. The universality of the chemical potential should be checked 
in the case of other noncritical strings using the two-matrix model, for example. 
Another interesting matrix model is the two-cut model, which corresponds to super 
Liouville theory. '^"^^ These problems arc left for future studies. 

In this paper, we have been attempting to determine the fundamental degrees 
of freedom in the nonperturbative formulation of string theory. It is suggested that 
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we should incorporate matrix, or at least something that has cndpoints, such as 
an open string. So far our investigation has been limited to the noncritical string 
theory case. However, it is conceivable that a similar situation exists for critical string 
theory. Should this issue be settled, it would be a great leap toward answering the 
ultimate question, What is string theory? 
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Appendix A: </)^-theory 

In this appendix, we consider one Hermitian matrix model with the partition 
function 



Zn = J d(j) exp(-Artry(<^)) , 



V{x) = lx'-y\ (A.l) 

to explain the convention used in the text. The critical value of the coupling constant 
in this theory is 

9. = ^- (A.2) 
2-34 

We define the square root part of the resolvent R{z) in the null instanton sector by 

(Wo{z)y^\[{V'{z))' + fo{z)) , (A.3) 

where fo{z) is a polynomial of degree 1. For g < g^, (Wo(x))^ has three zeros related 
as < X(2) < a^(o)- The quantity (Wo(x))^ is negative on the interval (iC(i), 3^(2)), 
where the eigenvalues are distributed continuously in the large- TV limit. At X(^q^, the 
first derivative of (Wo(x))^ vanishes . For g — g^, X(2) ^'^id X(o) coincide at 

x. = ii^i, (A.4, 
6 5* 
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while approaches 

V3-1 1 

^(i)(5*) = 5 — — ■ (A.5) 

We would like to determine (some ratios of) the constants P, 7 and k in the 
double scaling limit (4.6), so that the string equation and the sphere contribution to 
the free energy take the forms 

FWU„ = -l±*i. (A.6) 

After integrating over the angular variables of (p, Zq can be written 

Zo = Y[ dfi{xj) (/^({xj}j=i,...,7v))^ , 

d^i{x) = dxe-^^^''\ (A.7) 
where ^({a;j}j=i,- -,iv) is the Vandermonde determinant, 

^i{xj}j=i, -,N) = Hixi - xj) . (A.8) 

Prom the normalization of Pn{x), it can be written as 

^i{xj}j=i, -,N) = ^ det Pi_i(x,) . (A.9) 
By using this and the orthogonality of {-fn}„gz+u{o}' ^0 becomes 

N 

Zo = N\l[ hn-l 
n=l 

N-1 

n=l 

= exp ( In A^! + ^ (AT - n) Inrn I , (A.IO) 



' n 



n=l 
/ N-l 



\ n=l / 

where r„ is defined by 

rn^^ (nGZ+). (A.ll) 
By subtracting InA^!, the free energy F = lnZQ takes the form 

F= ^(A-n)lnrn. (A.12) 

n=l 
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In order to compute F, we need to know {r„}. From the equations 
n/in+1 = / dn{x) Pn-\{x) , 

= /rfa:^(e-^^(^)(P„(x))^) , (A.13) 

the recursion relations for {vn} and {^n} can be derived: 

ft 

rn{^- g{Sn + Sn+l)) = , 

sl-^Sn + rn + rn+i = 0. (A.14) 

The existence of the weak couphng limit g —>■ selects one of the solutions of the 
second equation, namely 

= l^^^HMr^+i^ . (A.15) 

When r„ approaches a continuous function of ^ = = 0(iV°) in the large-iV limit, 
the first equation in (A.14) implies that 

In the double scaling limit (4.6), we scale ^ and r(^) as 

KO = Y^e-«-«(^). (A.17) 

By applying Eqs. (A.15) and (A. 16) to the first recursion equation in Eq. (A.14), we 
obtain the Painlevc (I) equation governing the singular dependence of the specific 
heat M on t in the limit (4.6), (A.17): 

Also, the free energy in Eq. (A. 12) becomes 

1 AB'^T 

F{t) = ——^ ds{t-s)u{s). (A.19) 
9s 1^ Jt 

Prom Eq. (A. 18), the sphere contribution to u{t) reads 

«WUo = ;|/JVt. (A.20) 
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By inserting this into Eq. (A. 19), the universal part of the sphere contribution of 
F{t) is found to be 

By requiring the string equation (A. 18) to take the form given in Eq. (A.6), we 
obtain 

/3 ~ 3 ' 

We note that the same combinations of k, (3 and 7 appear in Eq. (A. 21). By using 
the values in Eq. (A. 22), it is easy to see that the sphere contribution necessarily 
takes the form given in Eq. (A.6). The following are useful expressions derived from 
Eq. (A.22): 

2 1 

«2 = -^/3i. (A.23) 

Next, we would like to take so that (Wo(C))^ approaches the conventional 
wo{C,) in Eq. (4.29) in the double scahng limit. If we choose 



13 3 , . /77, 9 

(Wo(x))^ approaches 



^ = - (2 + V3) = - {g^x^f , (A.24) 



{Wo{x)f = a^xZw {woiOf + Oia") , (A.25) 

where Z^y is given by 

Zw = l'x = eiH9*x*? ■ (A.26) 
Appendix B: 0^-theory 

Here we treat the Hermitian one matrix model with a quartic potential that is 
invariant under the transformation (p —(f), 

V{x) = \x'-^-x\ (B.l) 
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The value of the critical coupling constant is 



.. = ^- (B.2) 



{Wo{x))'^ in 0^-theory has four zeros positioned symmetrically: ±a;(o), ia^(i) (0 < 
< X(o))- For the one-cut solution, a cut runs along (— x^i), Assuming the 

existence of the weak coupling limit, X(o) and are given by 



^(0) 



^1(2 + ^1312^) 



^(i) = y^ (1-125)- (B.3) 

For g — > gi^i,, X(o) and come to coincide at 

x^ = 2V2. (B.4) 

Z2-symmetry implies that the coefficients s„ appearing in the recursion relations for 
Pnix) vanish. Thus, the recursion relations for r„ are found as 

Tl 

= J'n (1 - g{rn-l + Tn + r„+i)) . (B.5) 

We would like to obtain the form of this equation in the double scaling limit (4.7), 
together with 

g^ = g^e--"^', 
r(C) = ^e-«-«(^), (B.6) 

where C ^ ^■ 

In contrast to the (/>^ potential, the (j)^ potential has two local maxima, which 
arc located symmetrically with respect to = 0. Thus, it is necessary in the double 
scaling limit to focus on both maxima simultaneously. Therefore, we fix (3, 7 and k 
in Eq. (4.7) and r in (B.6) so that the string equation and the sphere contribution 
to the free energy take the forms 

t = Wt))'-f 8|«(t), 
Repeating the calculation given in Appendix A leads to 
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This implies that 

«^ = ^/3t. (B.9) 
We also demand that the continuum limit of {Wq{x))'^ approaches {wq{Q)'^: 

2 



(Wo{x)) =a^Zw (woiOr 

Vt 

Zw = l'x(^y. (B.ll) 



(55o(C)f=(c + v^) (C-^) , (B.IO) 



where ^ 



^3, 

Such a requirement determines the ratio as 

4 = 4. (B.12) 

T 

Appendix C: Behavior of r and s 

The orthogonal polynomials satisfy Schwinger-Dyson-type (SD) equations, 

= y"dx^ (e-^^(")Pn(x)P„(x)) . (C.l) 

The above SD equations can be translated into a set of recursion relations for the 
coefficients r„ and s„ appearing in Eq. (3.17): 

fn{r, S) = -5: 

. (C.2) 

, 9n{r, s) = 

It is easy to see from their derivation that the above equations contain the quantities 
r„ and s„ in such a way that this set of equations is invariant under the substitutions 



When we take n and N to infinity, fixing their ratio as 



(C.3) 

^n+j ^ ^n—l—j 
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the above-mentioned fact lead us to introduce continuum functions that correspond 
to r„ and Sn 




(C.5) 



with a possible 0{jj^) correction. With these definitions, the large-A^ limit of (C.2) 
becomes consistent and well-defined up to O (;^)- This can be seen by considering 
the fact that a term invariant under (C.3), say Sn+j + Sn-i-j, yields 2s{jj) + 
and no O (-^) correction. Substituting the coefficients r„ and s„ for the functions 
r(^) and Eq. (C.5) can be rewritten as 

—siCi / 1 \ 

«. = KO + ^ + 0(^), (C.6) 

if they are differentiable. Equation (C.2) determines r(^) and s(^) at each order in 
the ;^-expansion. 

It is a quite standard procedure to determine r(^) and s(^) explicitly. The 
large- AT limit of (C.2) can be summarized in the simple form 

where the contour of the integral is taken around the origin. As we tune the coupling 
constant g in V{(j)) to the critical value g^,, r(^) and s(^) behave as 

r-Tc^ {g*-gO^ > 

In the following, we give a more detailed description of the behavior of r(^) and s(^) 
near the critical region by using the expression (C.7). 

Equation (C.7) yields directly the following relations for f^{r,s) and g^{r,s) as 
functions of r and s 
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At the critical point, on the other hand, takes its maximum vahie 1 under the 
constraint = 0. Hence the Jacobi determinant for and vanishes: 



d{r, s) 

Substituting (C.8) into (C.9) yields 



0. (C.9) 



and further 



i^i) ---(I*) =0, (CIO) 

where we have taken the negative branch for the square root. The constraint = 
yields 

^ = -f^. (C.12) 

* 87-9^ 

Thus, we obtain the following relation at the critical point and for ^ = 1: 

^ = V^c. (C.13) 

Incorporating the above relation, the critical behavior of r and s can be written in 
the form 

r = rc - ^are Vl - ^ + (C.14) 

s = Sc - ^aV^Vl - C + 43, (C.15) 

where Ag = and a is a certain constant that might depend on the details of 
the model. The extra factor i accompanied by a is chosen so that (3.49) takes a 
familiar form. 

Appendix D: Pn{x) in the oscillating region 

In this appendix, we evaluate Pn{x) in the oscillatory region, applying analysis 
similar to that of the WKB method. In order to examine the continuity condition of 
the WKB method in the present context, it is first necessary to know the location 
of the turning point (s) for a specified x. We recall that r(^) is a monotonically 
increasing function of ^. In 0^-theory, s(^) is given by 
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Thus, s(i^) is also an increasing function of ^. For a given ^, we define as the 

zeros of the square root of q{x, 

X(±)(0 = 5(6±2V^. (D.2) 

The interval is the range of x on which q{x, ^) becomes purely 

imaginary. It is straightforward to see that increases monotonically from to 
a:(2) = = 1) > and that X(_) decreases monotonically from to = 

= 0) < 0. Such behavior of reveals the following simple structure. If 
X lies inside the cut 3^(2)], where the eigenvalues are distributed continuously, 

there is a no(^) such that 

^ ^ (77) ' ^(+^ (^i") < < AT . (D.3) 

If X lies outside the cut [x^i), x^)], there is no such no(x), so that q{x, is real. 

To summarize, we have learned that we should distinguish the two cases in which 
(a) X lies outside the cut [a;(i), a^(2)]) and (b) x lies inside the cut 3^(2)]- In case 
(a), there are no turning points, and therefore we only need the expression (3.28) for 
Pn{x) and apply it to the estimation of det(x — (f))'^- In case (b), there is a unique 
turning point ^o(a^) = ""jy'^ for a given x, and we need to obtain an expression for 
Pn{x) for no(x) <n<N. We next focus on this issue. 

Let us recall the continuity formula of the WKB method. As explained above, 
( < ^o{x) is the classically forbidden region, and ^ < (o{x) is the allowed region. 
Thus, the continuity formula for two independent solutions reads 

-J= e.p (i £ « ^ a 11 W«')l< + , (D.4) 

_i_ (-1 £ „,r)<) « ^ (i £ WOK . f ) . (D.5) 

Here, «;(^) is the wave number, and typically takes the form ~ ^ ^r an energy 

E and potential L'"(^) in a quantum mechanical system. Note that there emerges a 
factor of 2 in front of the sine on the right-hand side of Eq. (D.4), while Eq. (D.5) does 
not possess such a factor. In the following, we show that we need a multiplicative 
factor of 2 when we continue the expression of Pn{x) (3.28) analytically into the 
oscillatory region. For simplicity, we consider the case s = 0, that is, that in which 
the potential V{x) consists only of polynomials of even power, but the following 
argument also applies to the case s 7^ 0. 
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For a given x, there is a certain n(o) at which the argument of the square root 
changes sign in the expression of P„ (3.28): 

- 4r„(,, = 0. (D.6) 

In terms of the continuum variables = and ^(q) = ~W'^ ^^^^ point ^(q) corresponds 
to the turning point in the WKB method. Up to this point (with — 4r > 0, ^ < 
^(0)), the leading-order behavior of P„ in (3.28) is given by 

Pn ~ exp (^N j^^ log + i ^/x^ - 4r^ j (D.7) 

~ exp (-^^r'(C(o))(-r/)3) , (D.8) 

where we have introduced the new variable rj = ^ — ^(q) ■ 

We postulate that near ^^q)) both for 77 > and r/ < 0, P„, or the normalized 
polynomial V'n = -k=Pn, behaves as 



- sin (^-iV??!) or e^^-'')^ (D.9) 
Now, the quantities tp^ satisfy the recursion relation 

Xtpn = -v/rvPl'i/'n+l + V^V'n-l- (D.IO) 

We solve this recursion relation in the continuum limit near the turning point ^(g) 
and show that the postulated forms given in (D.9) are the actual forms. 
First, from (D.9), the difference of ^ behaves as 

|;* = ^|'/'~'<"V (D.ii) 

Thus for ?7 ~ , the higher-order differentials are suppressed. Expanding the recur- 
sion relation (D.IO) both in and ~ 0, we obtain the equation 

7ri + ^^^V' = 0. (D.12) 

This is the same equation as that which appears in the case of WKB method. Dis- 
regarding the factor '^^l'^' (or setting it to 1) for the sake of solving the equation, 
we have 

+ ryV = O. (D.13) 
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The following Airy function ^{ij) satisfies the above equation: 



1 



'm Jc 

The asymptotic form of ^{rj) for rj — oo is 

1 



1 r 

exp ( r]t + -t' 



— exp 

rj4 



2/ ^^ 



(D.14) 



(D.15) 



This form matches (D.8) if we set ^^^^ to 1. In the other asymptotic region, ri oo, 
<P becomes the trigonometric function 



1 



#(?7) - 2— sin { --(?7)5 + - ) . 



(D.16) 



Prom this, we conclude that there should be factor of 2 in front of the sine function, 
as in (D.16), when we continue (D.8) beyond the turning point. Also, both of the 
above asymptotic forms are consistent with (D.9). 

Appendix E: Evaluation of (^det{x — 0)^) for x inside the cut and a 
simpler evaluation of f^*^' da; e"^^''^^'"^ 

In this appendix, we present a calculation that allows for the evaluation of 
(det(x — ^)^) when x lies inside the cut. We also present a simpler evaluation of the 
integral of e~^^^^^^^ over the cut. 

First, let us evaluate (det(a; — (f))^) inside the cut (3.38): 

N N-l-i 

<det(x-^)2)= j2 (^^(^))' n ^^-i- 

i=no{x) j=0 

The expression for {Pi{x))^ is given by (3.31) . The other factor, Ylf=Q~^ rN-j, is 
evaluated as 

N-l-i 



N-l-i 

n ^^-j 

j=0 



exp 



exp 



exp 



j=0 



N 



N 



n I 1 



1 



d^lnr(0 + C 



1 



d0^r{^) + — Inr(l) - Trrrlnr — : 



2A^ 



2A^ 



N 



xll + O 



(E.l) 
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Hence, each term in Eq. (3.38) becomes 



N N-l-i 
i=no{x) j=Q 



=no{x) 

X , exp 
(^) 

"(0) 



2NRe I d^ln;fc(°)(x, ^ 

^0 



AT I'dC'lnriC) 



X exp 



2iV / d^lnk^^\x, 0+N / c/^lnr(0 

(E.2) 



where we have used ^k^^'' {x, — ''^(0 when ^ is greater than Now, we claim 
the relation 



7rp{x) 



1 



1 + °!^ 



This is proved as follows. Because Pn{x) = (det(a; — (f))), we have 

(exp [tr In (x - ^)]) = (det(a; - (f))) = Pn{x) 
A;W(x, 1) r .1 



(E.3) 



q{x, 1) 
X ( l + O 



exp 



iV / ci^lnfc(°)(a;, ^ 

^0 



Differentiating both sides of this equation with respect to x, we obtain 



(E.4) 



tr 



N \x- 



1 



det(a; 



(det(x - </.)) d^dMk^'^Hx, + djn 



(E.5) 



q{x, J J 

In the large-A?^ limit, the average on the left-hand side factorizes. Hence, Eq. (E.4) 
reduces at leading-order in N to 



Here we recall that 



ltr-l^)^|.ea.ln.(0)(.,o(l + C.(^)). 



(E.6) 
(E.7) 
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Thus, the imaginary part of this relation gives 

-2T,ip{x) = /l-tT—4—-) - /^tr ^ 



'0 

from which we obtain 



N X + ie — (p / \N X — ie — 



Prom the monotonic behavior of a;(±)(0, we have seen that there is such an no{x) 
that 

q{x + ie,0 = qix-te,0, (o < ^ < ^ 
q{x ± ie, = ±^ 61 • < C < l] (E.IO) 

Inserting this into Eq. (E.9), we obtain the relation in question, 

p(x) = - C di , / (l + O \ . (E.ll) 

Using the above relation in Eq. (E.2) and introducing ^o(a^) = , we obtain Eq. 
(3.39). 

Next, we show that there is a much simpler way to compute the left-hand side 
of (3.51) , 

■^den 

directly, without using the explicit form of Vesix). 

The point here is that we can replace the interval of the integration appropriately 
with a relative error within 0{1/N). First, as we saw in §3.3, the integrand is 0{1) 
inside the cut, while it decays exponentially as ~ exp(— A/"/(x)). Therefore, we can 
replace (E.12) as 

/oo 
(E.13) 
-OO 

because the contribution from outside the cut gives only cxp{—Ncx) ~ 0(1/N), 

where c is a certain constant of 0{1). Then, using (3. 7), (3. 14) and (3.15), we obtain 

Zden = {N+ l)hN ~ NhN. (E.14) 



dxe-^^-s^'^^ (E.12) 
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Next, we note that from (3.7), is obtained as 

/oo 
PNixfe-^^^'^Ux. (E.15) 
-oo 

However, we can again replace the interval of this integration 'inversely' with 

hN= PNixfe-^^^'^Ux, (E.16) 

up to a relative error of 0(1/N), for the same reason as above. Now, using (3.31) 
with n = N and noting that 

/•I 

V;5J) (x) = 2 Re / d^ln k^°^ (x, - V{x) (E.17) 
Jo 

is constant inside the cut, as we have seen in §2.1, we obtain 

hN = 2J^e-^^eH(^m) / , / ^ dx. (E.18) 

Finally, let us compute the remaining integral, 

/= / — -dx. (E.19) 

As mentioned in §§3.2 and 3.3, the cut is the region for which q{x,()^ < 0, and, 
in particular, q{x,l) provides the largest cut, < x < X(2)- Therefore, q{x,l) = 
\J{x- X(]^){x - X(2)), and 

'"^(2) 1 

dx = TT. (E.20) 

Hi) y(a;-X(i))(x(2) -x) 
Substituting this into (E.18), we obtain 

/iAr = 27rV^e-^^eff (^(2)) (E.21) 

and 

Zden = 2NT,^/HX)e-''^^S'(H2))^ (E.22) 
which of course agrees with the result in (3.51). 

Appendix F: Chemical potential for 0^ and 0^ potentials 

In this appendix we present explicit calculations of the chemical potential for 
the and potential cases. These two cases yield the same result as the other. 
This reveals the universality of the chemical potential. 
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The (^^-potential 

In this case, we use the potential V{x) = — |a;^. The resolvent is written 



1 



(F.l) 



The chemical potential /j, is given by (2.19). Hence, using (3.42) and (3.43), we 
obtain 

„(l-inst) 



,(0-inst) 



^ \{l - 2gs{l)f - Ag^r{l) 

((1 - 25^(1)) + V(1-25K1))2-45M1) 



■ exp 



ivi;(;)(x(o)) + ivv;S?)(x(2): 

(F.2) 



8Vr(l)((l-25s(l))2- 45^(1)) 

Here, we have used the saddle point method around x^q) = \~^{^) for the integration 
in the numerator of (2.19). 

We introduce t so that the free energy can be expressed in the universal form 



2 (gi-g 



2X§ 



gl 



4 5 
— t2. 

15 



15 V 3 

Then, the chemical potential ji can be written in terms of t as 



;7(1— inst) 



^(0-inst) g . 33/401:^1 



e 5 . 



(F.3) 



(F.4) 



The (/>^-potential 

In this case, the potential is given by V{x) 
1 Urn 



^x^ — ^x^. The resolvent is written 



R(x) = - ^V'{x) + (c/x^ + 2gr{l) - l) - 4r(l) 



(F.5) 



Using (3.42), (3.43) and (2.19), we obtain 



,7(l-inst) 
y(0— inst) 



1 



nN \ {l-2gr{l)){l~6gr{l)) 
X sVKl)(l-6.r(l)) [-^KS?^(-(o)) + , (F.6) 
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where a:(o) = y ^ " 2^(1). 

We introduce t so that the free energy can be expressed in the universal form 



r2^ f9c-9\^ , 4 ^5 



F = -N'—{^^^-^\ =-2—t2. (F.7) 
15 V 5c y 15 ^ ' 

Here, the factor 2 reflects the fact that there are two critical points at x = and 
X = —Xi,, because the potential is even. Then, we obtain the following value for the 
chemical potential: 



10 
11 
12 
13 

14 

15 
16 
17 
18 
19^ 
20 
21 
22 



^ = \^^ = 1 ^e-¥*^ (F.8) 
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